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ABSTRACT 

We report the Fermi-LKT discovery of high-energy (MeV/GeV) 7-ray emission positionally consistent with 
the center of the radio galaxy M87, at a source significance of over IOct in ten-months of all-sky survey data. 
Following the detections of Cen A and Per A, this makes M87 the third radio galaxy seen with the LAT. 
The faint point-like 7-ray source has a >100 MeV flux of 2.45 (±0.63) x 10~ 8 ph cm" 2 s" 1 (photon index = 
2.26 ± 0.13) with no significant variability detected within the LAT observation. This flux is comparable with 
the previous EGRET upper limit (< 2.18 x 10~ 8 ph cm" 2 s" 1 , 2a), thus there is no evidence for a significant 
Me V/GeV flare on decade timescales. Contemporaneous Chandra and VLBA data indicate low activity in the 
unresolved X-ray and radio core relative to previous observations, suggesting M87 is in a quiescent overall 
level over the first year of Fermi-LAT observations. The LAT 7-ray spectrum is modeled as synchrotron self- 
Compton (SSC) emission from the electron population producing the radio-to-X-ray emission in the core. The 
resultant SSC spectrum extrapolates smoothly from the LAT band to the historical-minimum TeV emission. 
Alternative models for the core and possible contributions from the kiloparsec-scale jet in M87 are considered, 
and can not be excluded. 

Subject headings: galaxies: active — galaxies: individual (M87) — galaxies: jets — gamma rays: observations 
— radiation mechanisms: non-thermal 
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1. INTRODUCTION 

As on e of the near est radio galaxies to us (D = 16 Mpc is 
adopted: lTonry|[l99"ll) . M87 is amongst the best-studied of its 
source class. It is perhap s best known for its exceptionally 
bright arcsecond-scale jet (Curtis 1918), well-imaged at radio 
through X-ray frequencies at increasingly improved sensitiv- 
ity and resolution over the decades (e.g. , Biretta et al l 1 199 lb 
iSparks et alj [19961 Marshall et alj 120021: iPerlman & Wilson! 
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120051) . Near its central ~ (3-6) x 10 9 solar mass supermas- 
sive b lack hole dMacchetto et all 1 19971: Gebhard t~& Thomasl 
2009), the jet base has been imaged down t o ~0.01 pc resolu- 
tion (~ 15-30x the Schwarzschi ld radius. Ijunor et al.lfl999b 
iLy et alJr200HlAcciari et"aT1l2009l) . 

At the highest energies, M87 is regularly detected by HESS, 
MAGIC, and VERITAS with variable T eV emission on 
timescales of years and flaring in a few days dAharonian et al.l 
I2006t lAlbert et alj[2008t lAcciari et al.ll2008 , 2009j). The sen- 
sitivity of these Cherenkov telescopes have also enabled the 
detection of another w ell-known nearby radio galaxy Cen A 
(lAharonian et al.ll2009l) . Without comparable imaging resolu- 
tion to the lower energy studies however, variability and spec- 
tral modeling are necessary to infer the production site of the 
TeV 7-rays and to deduce the source physical parameters. 

At high energy 7-rays (~20 MeV - 100 GeV), we are 
similarly poised for new radio galaxy discoveries with the 
Large Area Telescope (LAT) aboa rd the recently laun ched 
Fermi Gamma-ray Space Telescope dAtwood et al.ll2009l) . In- 
deed, we report here the detection of a faint, point-like 7- 
ray source positionally coincident with M87 using the Fermi- 
LAT. A fter the confirmation o f the EGRET disco very of 
Cen A (ISreekumar et all Il999h lAbdo et al l l2009d and in 
preparation), and t he recent detection of Per A/NGC 1275 
( Abdo et al. 2009b), this is the third radio galaxy successfully 
detected by Fermi. Unlike the known variable TeV source, 
there is so far no evidence for variability of the MeV/GeV 
emission in M87. An origin of the LAT emission from the un- 
resolved parsec scale jet (hereafter, denoted as the 'nucleus' 
or 'core') observed contemporaneously with Chandra and the 
VLBA 57 is discussed. Potential contributions from the larger- 
scale ( £ 0.1 — 1 kpc) jet to the unresolved 7-ray source are 
also briefly considered. Section 2 contains the details of the 
LAT observations, including a description of the Chandra and 
VLBA data utilized, with the discussion of these results in 
section 3. 



2. OBSERVATIONS 

The Fermi-LAT is a pair creation telescope which covers 
the en ergy range from ^20 MeV to >300 GeV dAtwood et al.1 
2009). It operates primarily in an 'all-sky survey' mode, scan- 
ning the entire sky approximately every three hours. The ini- 
tial LAT detection of M87 resulted from nominal processing 
of 6-months of all-sky survey data, as was applied to the ini- 
tial 3 -month dataset described i n I Abdo et al. (2009a), with a 
test statistic (M attox et al.ll 19961) . TS ~ 60. Including here an 
additional 4-months of data, the TS increased to 108.5, which 
is equivalent to a source significance ~ ^JTS = 10.4cr. The 
resultant 10-month dataset (Aug. 4, 2008 - May 31, 2009) 
corresponds to mission elapsed times (MET) 239557418 to 
265420800. Our analys is followed standar d selections of 
"Diffuse" class events ( At wood et al.1 [2009) with energies 
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FIG. 1 .— VLA A=90cm radio image from lOwen etaTU SOOO) with the LAT 
7-ray localization error circles indicated: rgs% = 5.2' and = 3.2' (statisti- 
cal only; see §|2). The M87 core is the faint feature near the center of the few 
kpc-scale double-lobed radio structure (in white). At the adopted distance 
D= 16Mpc, 1' =4.7kpc. 

E >200 MeV, a zenith angle cut of <105°, and a rocking 
angle cut of 43° applied in order to avoid Earth albedo 7- 
rays. Fermi Science tools 58 version v9rl0 and instrumental 
response functions (IRFs) version P6_V3_DIFFUSE were 
used for the analysis. 

A localization analysis with GTFINDSRC resulted in a 
best-fit position, RA = 187°.722, Dec. = 12°.404 (J2000.0 
equinox), with a 95% confidence error radius, 7-95% = 0°.086 = 
5.2' (statistical only; = 3.2'). To account for possible con- 
tamination from nearby sources, the model included all point 
sources detected at > 5er in an internal LAT 9-month source 
list within a region of interest (ROI) of r=15° centered on the 
7-ray position. Galactic diffuse e mission was modeled us- 
ing GALPROP ( Strong et aTll2004l) . updated to include recent 
gas maps and a more accurate decomposition into Galactocen- 
tric rings (galdef ID 5 4_5 9varh7S). An additional isotropic 
diffuse component modeled as a power-law was included. 
Figure Q] shows the re sultant 7-ray sourc e localization on a 
VLA radio image from lOwen et al.l (120001) . The 7-ray source 
is positionally coincident with the known rad io position of th e 
M87 core (RA = 187°.706, Dec. = 12°.391;|FejEitiD[2()04), 
with an offset (0°.020 = 1.2') that is a small fraction of the lo- 
calization circle. Currently, the best estimate o f the systematic 
uncertainty in 7-95% is 2.4' (lAbdo et al.ll2009ah . which should 
be added in quadrature to the determined statistical one. 

Spectral analysis was performed utilizing an unbinned like- 
lihood fit of the >200 MeV data with a power-law (dN/dE oc 
E~ r ) implemented in the GTLIKE tool. This resulted in 
F(>100 MeV) = 2.45 (±0.63) x 10" 8 ph cm" 2 s" 1 with a 
photon index, V = 2.26 ± 0.13; errors are statistical only. The 
flux was extrapolated down to 100 MeV to facilitate compar- 
ison with the previous EGRET non-detection of < 2.18 x 
10~ 8 ph cm" 2 s" 1 (2a ) from observations spanning the 1990's 
dReimer eT al. 2003). Thus, there is no apparent changes in 
the flux (i.e., a rise) in the decade since the EGRET obser- 
vations. Systematic errors of (+0.17/-0.15) x 10~ 8 ph cm" 2 

58 http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/Cicerone/ 
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FIG. 2. — The observed LAT spectrum (red circles) with representative 
TeV measurements of M87 in a low state from the 2004 observing season 
(black triangles) and du ring a high state in 2005 (blue squares), both by HESS 
lAharon ian et al. 2006). The lines indicate 1<t bounds on the power law fit to 
the LAT data as well as its extrapolation into TeV energies. 

s" 1 on the flux and +0.04/-0.1 1 on the index were derived by 
bracketing the energy-dependent ROI of the IRFs to values 
of 10%, 5%, and 20% above and below their nominal values 
at log(.E[MeV]) = 2, 2.75, and 4, respectively. The spectrum 
extends to just over 30 GeV where the highest energy photon 
is detected within the 95% containment. The LAT spectral 
data points presented in Figure [2] were generated by perform- 
ing a subsequent likelihood analysis in five equal logarithmi- 
cally spaced energy bins from 0.2-31.5 GeV. The 1 a bounds 
on the spectrum, obtained from the full >200 MeV unbinned 
likelihood fit, were extended to higher energies for compari- 
son with previous TeV measurements (see section 3). 

Lightcurves were produced in 10-day (FigureQ and 28-day 
(not shown) bins over the 10-month LAT dataset. Consider- 
ing the limited statistics, it was necessary to fix the photon 
index to the (average) fitted value in order to usefully gauge 
variability in the flux. Considering only statistical errors of all 
the binned data points with TS > 1 (lcr), a \ 2 test against the 
weighted mean fluxes of the 10-day and 28-day lightcurves re- 
sulted in probabilities, P(x 2 , v) = 22% and 70%, respectively, 
indicating plausible fits to the tested hypothesis. We conclude 
that there is no evidence for variability over the period of ob- 
servations. 

A radial profile of the 7-ray source counts (not shown) was 
extracted for the total energy range (>200 MeV). The profile 
is consistent with that of a point source simulated at ener- 
gies 0.2-200 GeV using the fitted spectral parameters above 
with a reduced \ 2 = 1 04 for 20 degrees of freedom. The total 
^0°.2 extent of the IP's kpc-scale radio lobes of M87 (Fig- 
ure [T] [Owen et a l. 2000) is comparable to the LAT angular 
resolution, 6» 68 ~ 0°.8 Ejg* dAtwood et al.ll2009t) . Therefore, 
from the presently available data, we can not disentangle (or 
exclude) a possible contribution of the extended radio features 
to the total 7-ray flux. 
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FIG. 3. — Lightcurve in 10-day bins obtained with the fitted photon index 
(r=2.26) fixed. The average flux is indicated with the solid horizontal line 
and the dotted lines are ±lcr about the average. Data points with TS < 1 (i.e., 
1<t) are shown as upper limits. 

To gauge the X-ray activity of M87 over the duration of 
the LAT observations, we analyzed five new 5 ksec Chandra 
ACIS-S images obtained in ^6 week intervals between Nov. 
2008 and May 2009 (PI: D. E. Hams). The X-ray core fluxes 
(0.5-7 keV) in these monitoring observations, (1.2 — 1.6) X 
1 0~ 12 erg cm" 2 s" 1 ( ~ 0.4-0.6 keV s" 1 in the units of Fig. 9 
of lHarris et al.ll2009h . are at the low end of the observed range 
over the last ^7 years. Additionally, the fractional variability 
is small (a7<flux>~ 0.1), indicating low X-ray activity in the 
core over the LAT observing period. 

At milli-arcsecond (mas) resolution in the radio band, M87 
has been monitored wi th the VLBA at 15 GHz since 1995 as 
part of the 2c m survey dKellermann et al.ll2004l) and MOJAVE 
dListeretalJ. 12009) programs . These data were re-imaged 
uniformly at 0.6 mas x 1.3 mas (position ang le=-ll°) reso- 
lution to match the additional map presented in K ovalev et alj 
d2007l) resulting in 23 total measurements of the unresolved 
core flux up to the latest observation on Jan. 7, 2009 (one of 
the Chandra exposures described above was obtained on the 
same day). This observation is the only one overlapping with 
the LAT dataset and the peak flux of 1.05 Jy/beam is consis- 
tent with the average over all the measurements (1.11 ±0.16 
Jy/beam). An indication of the sensitivity of these data to 
detecting flaring core emission is that the high flux state ob- 
served in the detailed 43 GHz VLBA moni toring at the time 
of TeV flaring in 2008 dAcciari et al J 120091 see § O is visi- 
ble in the 15 GHz data as a single high point on May 1, 2008 
(1.45 Jy/beam) 60 . This only suggests a period of low activity 
in the radio core (as reflected in the X-ray data), as the single 

59 See: http://www.physics.purdue.edu/MOJAVE/ 

60 Conversely, during the previous TeV flaring period spanning March - 
May 2005 I Aharonian et al. 2006), no comparable flare was visible in the 
VLBA 15 GHz core: 1.02 Jy/beam peak on Apr. 21st and 0.98 Jy/beam on 
Nov. 7th. Instead, the flaring X-ray/optical/radio knot HST-1 (60 pc from 
the core, projected) was observed to peak in early 2005 (Harris et al. 2009), 
suggesting an association with the variable TeV source I Cheu nget alJ2007r) . 



radio flux may not be representative of the entire 10-month 
LAT viewing period. 

3. DISCUSSION 

In blazars, it is commonly believed that 7-rays are pro- 
duced in compact emission regions moving with relativis- 
tic bulk velocities in or near the parsec scale core in order 
to explain the observed rapi d variability and to avoid catas- 
trophic pair-production (e.g., Dondi & Ghisellini 1995). Con- 
sequently , it is natural to exte nd this supposition to radio 
galaxies (Chia berge et al.ll2001l) which are believed to have 
jets oriented at systematically larger angles to our line of sight, 
thus constituting the parent population of blazars. Indeed, in 
the case of M87, a significant months timescale rise in the flux 
of the sub-pc scale radio core was discovered with the VLBA 
(at 43 GHz) during a period in e arly 2008 when fe w day 
timescale TeV flaring was detected dAcciari et al.l l2009). Dur- 
ing this period of increased activity, a Chandra measurement 
of the sub-arcsecond scale X-ray nucleus also ind icated a rel- 
ativel y higher flux than seen in past observations ( Harris et al. 
2009), thus signaling a common origin for the flaring emis- 
sions in the M87 nucleus. Therefore, during periods of lower 
7-ray activity, the radio/X-ray core can also be considered a 
dominant source of the unresolved higher-energy emission, 
and we discuss this in the context of the LAT MeV/GeV de- 
tection. 

In Figure 12 the LAT spectrum of M87 is plotted along 
with representative in tegrated TeV spectra from HESS 
d Aharo nian et al. 2006). The TeV measurements cover peri- 
ods when M87 was in its historical-mini mum (in 2004), an d 
during a high state (in 2005 - cf., Fig. 3 in lAcciari et al.l2 008). 
Although the formal difference in the fitted photon indices of 
the TeV data at high and low states is not statistically sig- 
nificant (r = 2.22 ±0.15 and 2.62 ±0.35, respectively), the 
LAT MeV/GeV spectrum (V = 2.3) connects smoothly with 
the low-state TeV spectrum. Taken together with the X-ray 
and radio measurements obtained during the LAT observa- 
tion (§ [2j, we view this as an indication that M87 is in an 
overall low 7-ray activity state during the considered period. 
In fact, no significant TeV flaring was detected in a prelimi- 
nary analysis of 18 hrs of cont emporane ous VERITAS obser- 
vations from Jan. - Apr. 2009 dHuill2009l) . 

M87 is the faintest 7-ray radio galaxy detected so far by 
the LAT with a MOO MeV flux (~ 2.5 x 10" 8 p h cm" 2 s" 1 ) 
about a n order of ma gnitude lower than i n Cen A ( Ab do et alj 
12009c!) and Per A dAbdo et all 12009b): the corresponding 
> 100 MeV luminosity, 4.9 x 10 41 erg s" 1 , is 4x greater than 
that of Cen A, but >200x smaller than in Per A. There is 
no evidence of intra-year or decade-timescale MeV/GeV vari- 
ability in M87 (§12), in contrast to the ^,7x and — 1.6x 
larger observed L AT fluxes than the p revious EG RET ones in 
the cas es of Per A dAbdo et al.l2009bl) and Cen A dAbdo et alj 
I2009cl) . respectively. The 7-ray photon index of M87 in the 
LAT band is similar to that of Per A {V = 2.3 and 2.2, respec- 
tively), while bei ng smaller than observed in Cen A (T = 2.9; 
Abdo et al. 2009c). These sources are low-power (FRI) radio 
galaxies, and have broad low-energy synchrotron and high- 
energy inverse Compton (IC) components in their spectral en- 
ergy distributions (SEDs) peaking roughly in the infrared and 
7-ray bands, respectively. Low-energy peaked BL Lac objects 
have similar shape d SEDs, with appr oximately equal apparent 
luminosities (e.g., iKubo et aT]|1998l) . As FRI radio galaxies 
are believed to c onstitute the parent popu lation of BL Lacs in 
unified schemes (Urry & Padovani 1995), the overall similar- 
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FIG. 4.— SED of M87 with the LAT spectrum and the Jan. 7, 2009 MO- 
JAVE VLBA 15 GHz and Chandra X-ray measurements of the core indi- 
cated in red. The non-simultaneous 2004 TeV spectrum described in Figure[2] 
and Swift/BXT hard X-ray limits (§[5J of the integrated emission are shown 
in lig ht brown. Historic al measurements of the core from VLA 1.5, 5, 15 
GHz iBiretta etaU1991l) . IRAM 89 GHz (Despringre et al.l 19961) . SMA 230 
GHz han et alJ|2008D . Svitzer 70, 24 u m IShi et alj|2007l) Gemini 10.8 /im 
IPerlman et alj|2001l). HST opti cal/UV ISparks et alJll996D . and Chandra 1 
keV from Marshall et al. 12002, hidden behind the new measurements) are 
plotted as black circles. The VLBA 15 GHz flux is systematically lower 
than the historical arcsec-resolution radio to infrared measurements due t o 
the presence of intermediate scale emission (see e.g., [Kovalev et afl l2007ji 
The blue line shows the one-zone SSC model fit for the core described in §[5] 

ity of their SEDs is not surprising. 

We construct a SED for M87 (Figure gjl using the the LAT 
7-ray spectrum and the overlapping Jan. 7, 2009 Chandra and 
VLBA measurements of the core. Also plotted are histori- 
cal radio to X-r ay fluxes of the core (see ISparks et al.l[l996l: 
iTan et all 12008) measured at the highest resolutions at the 
respective frequencies. The core is known to be variable, 
with factors of ~ 2 changes on months timescales common in 
the op tical and X-ray bands (Perlman et al. 2003; Har ris et alj 
2009). To help constrain the overall SED at frequencies be- 
tween the X-ray and LAT measurements, we determined in- 
tegr ated 3a upper lim its in three hard X-ray bands (follow- 
ing, Aiello ^etal] 120081) based on the Swift/BAT dataset in 
Aiello et all d2009l) 7 including about another additional year 
of exposure (i.e., ~4 years total from Mar. 2005 - Jan. 2009). 

The broad-band SED is fit with a homo geneous one-zone 
synchrotron self-Compton (SSC) jet model dFinke et alJ2008l) 
assuming an angle to the line of sight, 8= 10°, and bulk 
Lorentz factor, T\, = 2.3 (Doppler factor, S = 3.9), consistent 
with observations of appare nt motions of ^ 0.4c (Tb > 1.1) 
in the parsec-scale radio jet (Lv et al. 2007). A broken power- 
law electron energy distribution N(j) oc 7~ p is assumed, and 
the indices, p\ = 1.6 for 7 = [1, 4 x 10 3 ] and p2 = 3.6 for 
7 = [4 x 10 , 10 7 ] are best guesses based on the available 
core measurements. The normalization at low energies is con- 
strained by the single contemporaneous VLBA 15 GHz flux 
which is measured with ~ 10 2 - 10 3 x better resolution than 
the adjacent points. The source radius, r = 1.4 x 10 16 cm = 
4.5 mpc is chosen to be consistent with t he best VLBA 43 
GH z map reso lution (r <7.8 mpc = 0. 1 mas, Uunoret all 1991 
ILv et al] 120071) and is of order the size implied b y the few 
day timescale TeV variability dAcciari et al.l 120091) . For the 
source size adopted, internal 7-7 absorption is avoided so 
that the LAT spectrum extends relatively smoothly into the 
TeV band, consistent with the historical-minimum flux de- 



tected by HESS ( Aharonian et al. 2006) and the prelimi nary 
upper limit of < 1 .9% Crab from VERITAS observations dHuil 
2009) contemporaneous with the LAT ones. 

In the SSC model, the magnetic field is B = 55 mG and 
assuming the proton energy density is 10 x greater than the 
electron energy density, the total jet power is Pj ~ 7.0 x 10 43 
erg s . The jet power is particle dominated, with only a 
small contribution from the magnetic field component (Pb ~ 
2 x 10 40 erg s -1 ). In comparison, the total kinetic power in the 
jet is ~Tew x 10 44 erg s" 1 as d etermined from the energetic s 
of the kpc-scale jet and lobes dBicknell & Begelmanlll996l) . 
and is consistent with the jet power available from accre- 
tion^ < 10 45 erg s" 1 dRevnolds et al.lll996t iDTMatteo et alj 
120031) . These power estimates are similar to those derived 
for B L Lacs from similarly mo deling their broad-band SEDs 
(e.g. JCelotfi & Ghisellinill2008l) . 

As applied to M87, such single-zone SSC emission mod- 
els also reproduce well the broad-ban d SEDs up to MeV/ GeV 
energie s in the radio galaxies Per A dAbdo et alj|2009b1) and 
Cen A dChiaberge et al.ll20 01). In this context, the observed 
MeV/GeV 7-ray fluxes of b lazars appear to be correlated with 
their compact radio cores dAbdo et alj|2009ct IKovalev" et al. 
2009), suggesting a common origin in the Doppler boosted 
emission in the sub-parsec scale jets. The fact that the 3 radio 
galaxies detected by the LAT so far have amongst the bright- 
est ( ^ 1 J y) unresolved radio co res, in line with these ex- 
pectations (Ghisellin i et al.l 120051) . lend evidence for a com- 
mon connection between the 7-ray and radio emitting zones 
in such jets. 

It should be emphasized that these observations are not 
simultaneous and particularly, the TeV emission is known 
to be variable on year timescales, so other emission com- 
ponents may contribute to the variable emission. There- 
fore, although not strictly required, more sophisticated mod- 
els over the single-zone one presented can reproduce or con- 
tribute to the observed emission. In particular, the beam- 
ing requirements in the one-zone SSC modeling of the three 
known 7-ray radio galaxies are systematically lower than re- 
quired in BL Lacs, s uggesting velocit y profiles in the flow 
(Chia berge et al.ll2001l). Such mode ls dGeorgano poulos et alj 
120051: Tavecchio & Ghisellini 2008) have in fact been used 
to fit the SED of M87 in addit ion to models ba sed on 
additional spatial structure (e.g.. iLenain et alJ 12007). Pro- 
tons, being inevitably accelerated if they co-exist with elec- 
trons in the emission regions, probably dominate energet- 
ically and dynamically the jets of powerful AGN (e.g., 
Celotti & Ghisellini 2008). Applying the synchrotron -proton 
blazar model dMiicke et alj|2003t iReimer et al] 120041) to the 
quiescent M87 data set yields reasonable agreement model 
fits that support a highly magnetized compact emission region 
with approximate equipartition between fields and particles 
and a total jet power comparable with the above estimates, 
where protons are accelerated up to ~ 10 9 GeV 

Outside of the pc-scale core, the well-known arcsecond- 
scale jet (e.g.. IBiretta et all I199U Marshall etaT] 120021: 
Perlman & Wilson 20051) is also a possible source of IC emis- 
sion. As both the LAT and TeV telescopes are unable to 
spatially resolve emission on such small scales, the expected 
spectral and temporal properties of the predicted emission 
must be examined. On the observed scales, the dominant 
seed photon source is the host galaxy starlight, and such 
an IC/starlight model applied to one of the brightest re- 
solved knots in the jet - knot A, ^1 kpc projected distance 
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from the core - results i n a spectrum peaking at TeV ener- 
gies dStawarz et al. 2005), thus producing a harder MeV/GeV 
spectrum than observed by the LAT. Even closer to the core 
(~60 pc, projected ) the superlum inal knot HST-1 (>4c-6c; 
iBiretta et alJTT999l: (Cheung et all [2007) is a more complex 
case. This knot is more compact than knot A, and its IC 
emission is expected to be further enhanced by the increased 
energy densities of the surrounding circumnuclear and galac- 
tic phot on fields, as well o f the comoving synchrotron ra- 
diation dStawarz et al. 2006). The radio/optic al/X-ray fluxes 
of HS T-1 have been dec lining since its giant flare peaked in 
2005 (lHarris et al.l l2009). with current X -ray fluxes compara- 
ble to its pre-flare levels in 2002. Considering the variable 
and compact nature of the source (with observed months dou- 
bling timescales implying r £ 228 mpc; cf., footnote l60l. the 
predicted IC spectrum has a complex temporal and spectral 
behavior. In the absence of detailed contemporaneous mea- 
surements, its possible role in the production of the LAT ob- 
served MeV/GeV emission is unclear. 

Continued LAT monitoring of M87 coordinated with multi- 
wavelength observations can extend the current study of 'qui- 
escent' emission to possible flaring, in order to further ad- 
dress the physics of the radiation zone. While the extragalac- 
tic 7-ray sky is do minated by blazars (Hartman et al] 1 19991 : 
lAbdo e t al. 2009c), this optimistically indicates an emerging 
population of 7-ray radio galaxies. Other examples, includ- 
ing the few possible associations w ith EGRET d etections like, 
NGC 625 1 ilMukheriee et alJl2002h and 3C 1 1 1 (ISguera et al] 
120051: iHartman et al.ll2008l) await confirmation with the LAT, 
and more radio galaxies are expected to be detected at lower 
fluxes. This holds great promise for systematic studies of rel- 
ativistic jets with a range of viewing geometries in the high 



energy 7-ray window opened up by the Fermi-LAtT. 
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